Constraints on Shift-Symmetric Scalar-Tensor Theories with a Vainshtein Mechanism 

from Bounds on the Time Variation of G 
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We show that the current bounds on the time variation of the Newton constant G can put 
severe constraints on many interesting scalar-tensor theories which possess a shift symmetry and 
a nonminimal matter-scalar coupling. This includes, in particular, Galileon-like models with a 
Vainshtein screening mechanism. We underline that this mechanism, if efficient to hide the effects 
of the scalar field at short distance and in the static approximation, can in general not alter the 
cosmological time evolution of the scalar field. This results in a locally measured time variation of 
G which is too large when the matter-scalar coupling is of order one. 
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Many theories in which gravity is modified with re- 
spect to general relativity (GR) contain, in addition to 
the metric, a scalar field which is coupled directly to mat- 
ter. Such scalar-tensor theories appear naturally in low 
energy limits of string theory and are also obtained from 
phenomenological brane- world constructions (such as the 
DGP model [l|). Some are also of current interest as able 
to produce an interesting cosmology via a large distance 
modification of gravity. In such theories, in contrast to 
GR, matter not only interacts via the helicity-2 graviton, 
but also via the exchange of the scalar field. In general, 
one faces the following dilemma: Either this field is cou- 
pled to matter with gravitational strength, as required 
to produce order one deviations from GR, but then the 
theory cannot pass local tests of gravity, or the coupling 
is very small, but then there are no significant effects 
of the scalar. A canonical example is the Brans-Dicke 
theory [2, Q and its extensions fill whose parameters are 
tightly constrained by the local tests of gravity and ob- 
servations of binary pulsars (see for instance [5| ) . 

A way out of this dilemma is provided by the Vain- 
shtein mechanism, first proposed in the context of mas- 
sive gravity [y, |7| (a proof was recently provided in [8|). 
Indeed, close to localized bodies, this allows to screen 
effects which lead to large deviations from GR at large 
distances. This mechanism was also shown to be present 
in the DGP brane model [1] as well as its decoupling limit 
[9j] . It was later generalized and shown to apply to a large 
class of scalar-tensor models, called in [l3| "k-moufiage" 
gravity models, with a non-linear kinetic self-interaction 
of a scalar field providing a self-screening of the scalar 
force a la Vainshtein (hence the name k-moufla ge) . This 
class contains in particular the Galileon model |ll| . and 
its covariantized versions ^12h15| . Many applications of 
the Galileon model and its extensions to the late-time 
acceleration, including minimally coupled |16l4l9| as well 



as non- minimally coupled models [20|-|22| , have been con- 
sidered, while various constraints coming from cosmol- 
ogy as well as from local observations have been studied 

In this paper, we point out that in spite of the fact that 
the Vainshtein screening indeed allows to pass most of 
the constraints coming from local observations by cutting 
off the spatial variation of the scalar field near massive 
bodies, the tests on the constancy of the Newton constant 
may easily rule out many models. Indeed, we show that 
in many shift-symmetric models, the evolution with time 
of the scalar field is (approximately) the same everywhere 
and it follows its cosmological behavior. If the scalar is 
directly coupled to matter, this induces a variation of 
the Newton constant G, which is tightly constrained by 
a number of observations (see, e.g., the review [2g|). The 
most stringent bounds come from binary-pulsar data [27[ 
and above all Lunar Laser Ranging experiments [2^ , the 
latter giving \G/G\ < 1.3 x 10~^^yr~^, or in terms of 
the Hubble value today, Hq, \G/G\ < 0.02iJo. As we 
will see below, the time variation of the scalar field is 
generically of order of the Hubble scale Hq (unless it is 
in the "cosmological" screening regime with a tiny energy 
scale M <C Hq). This, whenever the direct coupling of 
the scalar field to matter is of order of one, induces a too 
large variation of Newton's constant, \G/G\ ^ Hq. 

We consider the following general action. 
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where R is the Ricci scalar of the metric g^j/ 
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Cs = 

— {dipy is the standard kinetic term of a scalar field f 
(normalized to be dimensionless) , £nl describes some 
generic nonlinear self-interaction of tp, and the matter 
fields (collectively denoted as ipm) are minimally coupled 



to the physical metric g^„ = A'^{(p)g^„. Because of this 
couphng, gravity is modified at large distances through 
a scalar exchange, while GR is supposed to be restored 
at small distances thanks the Vainshtein screening effect 
made possible by the self- interaction £nl- This screen- 
ing occurs for rather generic nonlinear interaction terms 
-Cnl [10] and we do not specify a precise from of it. We 
however assume that £nl is shift symmetric, i.e., it does 
not change under the transformation ip —> ip + const. For 
example, we can choose £nl to be in the Galileon or k- 
essence families. For dimensional reasons Z^nl niust con- 
tain a mass scale M. We will assume that this is the only 
additional scale entering our action. In general, this scale 
is fixed by phenomenological requirements, e.g. to get 
present day acceleration of the Universe from Galileons, 
M should be of order of the Hubble scale Hq. 

Note that by a conformal transformation, action ([T]) 
can always be rewritten in a form where matter is only 
minimally coupled to one metric 5^1,, going to the so- 
called Jordan frame. Here we will work rather in the 
Einstein frame ([T]), with the understanding that our re- 
sults would apply to any theory whose action can be put 
in the form ([T|) by a suitable field redefinition. We will 
also not consider effects that can arise, even in theories 
with minimal coupling to gravity in the Einstein frame, 
due to the non-linear kinetic coupling of the graviton to 
the scalar field (called "kinetic braiding" in (l6{). 

The variation of Eq. ([T|) with respect to the metric g^^i, 
gives the (modified) Einstein equations. 
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where T^t , TJii, ' and T^™ are respectively the energy- 
momentum tensors for the standard scalar kinetic term, 
its nonlinear term, and the matter contribution. The 
equation of motion for the scalar field is 



V^(VV + J^l) 
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(3) 



where a{ip) = din (A) /dip, the nonlinear current J^^ is 
obtained by variation of £nl with respect to the gradient 
of the scalar field, J^^ = ~^^^nl/S(p^i^, and T^™) is 
the trace of the matter energy-momentum tensor in the 
Einstein frame. Note that the field equations ([3]) can 
always be written as the divergence of a current, because 
of shift symmetry. We also rescale the Planck mass so 
that A{ip) — 1 at present. We stress that the trace of the 
matter energy-momentum tensor defined in the Einstein 
frame, T/i™% differs by the factor A'^{ip) from the trace of 
the (conserved) Jordan-frame energy momentum tensor. 
However, as we will see below, the time variation of if is 
small (of order of the Hubble scale or less) , so that the 
change of (p with time can be neglected in the r.h.s. of ([3]) , 
giving only small corrections. For instance, ii A{ip) = e"^ , 
the approximation A'^{tf) sa 1 is valid for \ip\ ^ 1, i.e., 
for |At| < i?-i. 

The cosmological evolution of the scalar field, (^cosm (0 
can be easily read from ([3]), 



Three different cosmological regimes can be identified. 
In a regime where the cosmological energy density of the 
scalar field, p^, is subdominant compared to the mat- 
ter energy density, p^, Ptp ^ Pm, the scalar field equa- 
tion Q is decoupled from the metric equation ([2]). Then 
from the Einstein equations it follows that pm — 3MpH^, 
thus the r.h.s. of (jlj is ~ aH^. If the scalar field is 
away from the "cosmological" Vainshtein regime (i.e., 
when the nonlinear term in the l.h.s. of ([4]) is negli- 
gible), then from Eq. ^ one can see that a particular 
solution to dH) is |(/3cosm| ^ OiH. Notice that the general 
solution also contains a homogeneous decaying solution 
Coexp (^—3 J dm) with an arbitrary constant Co, how- 
ever unless this constant (or, equivalently, the initial con- 
dition) is fine tuned, the time variation of (p remains of 
order of aH. 

In the second regime matter is again dominant, p^ <C 
Pni7 but the scalar field is in the cosmological Vainshtein 
regime. Then Eq. ^ contains, in addition to H, also the 
"nonlinear" scale M. Therefore the solution to ^ gener- 
ically contains a combination of scales H and M. When 
this scale is small with respect to Hq, the time evolution 
of V'cosm may thus be suppressed, |(/9cosm| *^ ^o- This is 
the cosmological analog of the original Vainshtein mech- 
anism. 

The third regime is realized in the case when the scalar 
field is dominant, p^ 3> Pm, and in particular when the 
late-time acceleration of the Universe is driven by the 
scalar field. In this case, both the metric and the scalar 
field equations depend only on one dimensionful param- 
eter, M, which is of order of Hq. We thus conclude that 
the typical value of the present variation of the scalar 
field is the Hubble scale, |(^cosm| ~ -^o- 

Let us consider now the local effects caused by the 
conformal coupling of the scalar field. For a slow time 
evolution, (pcosm(^) can be written as the linear approxi- 
mation, 



V'cosiiHO — Vcosm{ta) + 'Pcosm(io) ^- 



(5) 



Note that ^ imposes the boundary value of the field far 
from localized sources. The solution to the full equation 
of motion ^ at any point of space-time (including the 
regions close to massive bodies, in particular, inside the 
Vainshtein radius) depends on both time and space coor- 
dinates. The key observation here is that thanks to the 
shift symmetry of the equation of motion, the PDE ^ 
allows separation of variables in the following way. 



p{t, r) = p{r) + (pcosm(io)i + <(5cosm(io) 



(6) 
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where r is the distance to the source. It is not difficult 
to see that the above ansatz (which has also been used 
in other contexts ^3|) "passes through" the full equation 
of motion ([3]), giving an ordinary differential equation of 
the second-order on p{r), with possible remnants from 
the time-dependence in a form of a constant, (/'cosm(^o)- 
The last two terms in the above ansatz give the boundary 
condition for the PDE imposed by the cosmological evo- 
lution ^ , provided we choose the radial dependent part 



of (O, ip{r), to vanish at infinity. Now, the ODE on ip{r) 
is of the second order, and supphed with two boundary 
conditions, (f{r — oo) = and (p'{r = 0) ~ (the last 
one comes from the regularity at the origin), which is 
in general sufficient to find a unique solution. Provided 
that this solution is non singular (which is in some cases 
a strong mathematical requirement), and assuming that 
([S]) is a good approximation for the time-dependent cos- 
mological evolution of (^, our ansatz gives a solution for 
all times to the field equation ([3]) for a spherical source 
centered at r = 0. The key point of this paper is that the 
time derivative of ip{t, r) is found to be independent of 
r, i.e., it is set by the cosmological evolution even inside 
the regions where the Vainshtein screening operates. 

It is worth mentioning that our ansatz (|6]) and bound- 
ary conditions are in fact selecting a particular class of so- 
lutions. Indeed, if formulated in terms of a Cauchy prob- 
lem, they implies ip{t = 0,r) = ip{r) and ip{t — 0,r) — 
V'cosm(io)- In principle, there is no reason not to choose 
some radially dependent initial velocity, ip{r) = C{r). In 
contrast to ([6]), such solutions, however, are not station- 
ary, and they should relax to the stationary one, assum- 
ing that the latter is stable. A numerical check of this 
relaxation, as required by the nature of the field equa- 
tions, goes however far beyond the scope of this paper. 

Let us consider a couple of illustrative examples. First, 
as the non-linear term, we take one of the Galileon La- 
grangians, namely. 



£nl = -M-^Dip{d^)^ 



(7) 
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It is not diffi- 
cosm in two different 



and the coupling to matter A{<^) — e 
cult to find that the evolution of (p, 
cosmological regimes is in accord to our general findings 
(we assume here a vanishing cosmological constant): (i) 
when p^p <^ pm and the nonlinear term is subdominant 
in ((H) then i^cosm = — 2H; (ii) when p^p ^ pm and the 
nonlinear term is dominant in Q then iPcosm — 2Af ^/3. 
The time dependent solution for the scalar field around 
a body of mass m is then given by ([6]) with. 



^'(r) = 



rM^ 




(8) 



Note that the last piece inside the parentheses is due 
to the cosmological time evolution of (p. As a second 
example we consider a scalar field Lagrangian with the 
the signs of the scalar field kinetic terms fiipped with 
respect to previous example, and the coupling to matter 



Aiip) 



This Lagrangian allows self-accelerating 



solution [3, with H^ = MV(3^/6) and ip = V&H, while 
the time-dependent stationary solution for the scalar field 
is still given by (jS]) with the same radial-dependent part 

The time-dependence of ip leads to a variation of the 
effective Newton constant with time. This can be seen 
by making the conformal transformation to the Jordan 
(physical) frame, with the metric g^i^. Generically, there 



are two effects, which enter the final result for the evo- 
lution of Newton's constant: the exchange of helicity-0 
modes and the rescaling of the coordinates via the con- 
formal transformation of the metric. It should be noted 
that in standard scalar-tensor theories (without screening 
mechanisms), these effects arc of the same order, so that 
they even can compensate each others (as in Barker's 
theory [BOj ) giving no change of G. In our case, how- 
ever, the exchange of ip is screened by the Vainshtein 
mechanism, so that only one effect — the stretching 
of coordinates — is important. As a result, the effec- 
tive Planck mass in the action gains a dependence on 
'Pcosm(i), Mp = A~''- {(pcosm) Mp . Thus, the observed 
Newton constant evolves with time as 



\G/G\ « 2a(pcosm(i)- 



(9) 



As we have seen before, depending on the regime, 
|</5cosm| '^ OiH (in the matter domination regime) or 
|</5cosm| ~ H (when the scalar field is dominant). There- 
fore, presently one has \G/G\ ^ o^Hq if the scalar field 
is subdominant and away from the cosmological screen- 
ing, and \G/G\ ~ aHg when the scalar field is dominant, 
in particular, when it drives the late-time acceleration of 
the universe. This applies in particular for a constant a, 
i.e., a conformal coupling g^^i, = e^°"^gp,^. 

The observational constraints from Lunar Laser Rang- 
ing give \G/G\ < 0.02iJo which is enough to rule out 
theories of the kind considered here with a scalar cou- 
pling to matter of the order of the gravitational one (i.e., 
a ~ 1). In order for a theory to explain the acceler- 
ated expansion of the universe at present days, and pass 
the constraints on the variation of G, one should assume 
a < 0.01. It is interesting to note that a similar con- 
straint on the matter-scalar coupling constant was ob- 
tained for the covariant Galileon theory from a combined 
analysis of supernovae, baryonic acoustic oscillations and 
cosmic microwave background |22| . 

Let us now briefly discuss the case of non shift- 
symmetric theories. There is a class of such theories 
which can be put in the form ([1]) by suitable field re- 
definitions, in which case our conclusions apply. When 
this is not the case, the situation must be carefully re- 
analyzed. Indeed, the fact that the ansatz (jG)) leads to 
a mere ODE to solve for ip{r) is a direct consequence of 
the shift symmetry and it might be that some screening 
of the time variation of G occurs when this symmetry is 
lost. For example, we may introduce a mass term in the 
action, rn?ip'^, which explicitly breaks the shift symme- 
try. If the mass is big enough (say, much bigger than the 
present Hubble scale, Hq), then the cosmological evolu- 
tion of ip is suppressed, because the scalar follows the 



minimum of the effective potential, tpn 



aHH/r 



However, such theories do not possess either an interest- 
ing self-accelerating scenario driven by the kinetic term. 
We also stress that even in the shift-symmetric case, it 
might be that the ODE obeyed by (p{r) does not possess 
regular solutions, or that it leads to a stationary solution 
^ which happens to be unstable as a solution of the 



PDE ([3]). If so, it may open a way out of our conclusions, 
necessitating the use of a more general ansatz than © to 
solve the field equations, which could in turn result in a 
Vainshtein suppression of the time variation of G in the 
solar system. 

It is also interesting to mention another possibility of 
avoiding any significant time evolution of G, namely, to 
violate our assumption of a conformal matter-scalar cou- 
pling. In particular, in the relativistic MOND theory, 
called TeVeS [31| , where the physical metric is related to 
the Einstein one in a disformal way, the time variation 
of the Newton constant is strictly zero [33 ■ This is a 
consequence of the different scalings of time and space 
coordinates with A{ifcosm) when one imposes that the 
physical metric g^j,u tends to the Minkowski metric at 
spatial infinity. This absence of time variation of G also 
applies to other theories with disformal coupling, in par- 
ticular to the improved relativistic MOND [33] , where the 
k-mouflage screening has been used to pass solar-system 
and binary-pulsar constraints. 

Let us finally underline that the DGP brane model, 
although equivalent to a scalar-tensor theory of the 
Galileon type in the UV (in particular in the decoupling 
limit 19| ) , is not fully described by such a theory at cosmo- 
logical scales (IR limit). Therefore our analysis does not 
apply to the DGP model, and this explains why Ref. jsj 
did not find any local time variation of G. 

To conclude, we have shown that a generic scalar- 
tensor theory with conformal coupling of a scalar field 
to matter, and with a shift symmetry of the scalar La- 



grangian, is tightly constrained by the bounds on time 
variation of the Newton constant. The models which fall 
into this category are not only standard (massless) Brans- 
Dicke-like theories, but also those featuring a Vainshtein 
screening mechanism due to the kinetic self-coupling, in 
particular non-minimally coupled Galileon models. We 
argued that the local time evolution of the scalar field 
is set by its cosmological evolution and not screened by 
the Vainshtein effect, which is only able to suppress the 
"fifth force" due to the exchange of helicity-0 degree of 
freedom. The time derivative of the scalar field is of 
order of the Hubble scale, unless the whole universe is 
in the screening regime ( "cosmological" Vainshtein effect 
with M -^ Hq). This induces a time evolution of New- 
ton's constant of the same order. This result applies for 
both the matter domination and the scalar field domi- 
nation cosmological regimes. It also does not depend on 
a particular form of the non-linear self-interaction term, 
provided it is shift-symmetric. The key point is that the 
time dependence in ([B]), which is crucial in ©, is the 
same irrespectively of the precise structure of >Cnl • The 
evolution of the Newton constant is, however, tightly con- 
strained by observations, \G/G\ < 0.02_ffo, therefore the 
conformal coupling on such theories is constrained too. 
For example, if the non-minimal coupling is of the form 
ex-p{2a(p), then \a\ should be less than 10^^. 
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